Molecular nanomagnets as quantum simulators 
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We start by considering different kinds of one- 
dimensional, translationally invariant Hamiltonians TL 
that can be mapped onto a model of 1/2 spins s.^, with 
nearest-neighbour (NN) interactions and Hamiltonian H. 
Our aim is to simulate the time evolution associated with 
H ("target evolution" U{t)) by means of the proposed 
hardware. The first step is to approximate U{t) by the 
Trotter-Suzuki formula (h — 1): 
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where H = H^^L + H'^^\ The contributions H^^^ 
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(Hcvcn) include all the two-spin terms 
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Quantum simulators are controllable systems that can be used to simulate other quantum systems. 
Here we focus on the dynamics of a chain of molecular qubits with interposed antiferromagnetic 
dimers. We theoretically show that its dynamics can be controlled by means of uniform magnetic 
pulses and used to mimic the evolution of other quantum systems, including fermionic ones. We 
propose two proof-of-principle experiments, based on the simulation of the Ising model in transverse 
field and of the quantum tunneling of the magnetization in a spin-1 system. 

PACS numbers: 03.67.Ac 03.67.Lx 75.50.Xx 

The simulation of quantum systems by a classical com- 
puter is intrinsically inefficient, because the required 
number of bits grows exponentially with the system size. 
This makes many important problems in physics and 
chemistry intractable. Such limitation might be over- 
come by quantum simulators (QSs), whose dynamics can 
be controlled so as to mimic the evolution of the target 
system [l|. 

The implementations of quantum simulators so far pro- 
posed essentially fall into one of two categories 0] ■ In the 
first one, large and globally addressable systems are used 
to analogically simulate specific target Hamiltonians. In 
the second one, QSs consist of few, individually address- 
able qubits, and the time evolution of any target system 
can be discretized into a sequence of logical gates. Here 
we propose a hybrid approach, where the simulation of 
different kinds of translationally invariant models is per- 
formed by exploiting chains of molecular nanomagnets, 
manipulated by means of spatially homogeneous mag- 
netic fields. As in analog QSs, our hardware consists of 
a potentially large array of qubits, with a geometry re- 
flecting that of the target system, and no local control 
required. As in digital QSs, the manipulation of the QS 
state is here achieved by suitable sequences of quantum 
gates, performed in parallel on the whole array. This 
permits to simulate a large class of models for each given 
geometry. 

Our proposal exploits two classes of molecular nanomag- 
nets [Sj that play two distinct roles: (effective) S — 1/2 
spins are used for encoding the qubits; these are con- 
nected through antiferromagnetic systems with nonmag- 
netic ground state (3 = 0), that can be controUably ex- 
cited so as to effectively switch the coupling between the 
qubits. The capability of engineering complex structures 
consisting of weakly coupled and monodispersed systems 
has been recently demonstrated in molecular magnetism 
The wide tuneability of both the intra- and inter- 
molecular interactions, combined with the possibility of 
coherently driving the spin dynamics [RQ, makes these 
system suitable for both spintronics [7| and quantum- 
information [§-[1^ applications. In addition nanomag- 
nets can be g rafted onto surfaces without altering their 
properties [l3|. 
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(^cvcn(s2fc, S2fe+i)), whilc H^^^ includcs all the single-spin 

terms. Since H^^'^, Hcven, and Hi generally don't com- 
mute, Eq. (m is only exact in the limit r = t/n -> 0. 
Each of the three terms in parentheses in Eq. ([T|) can be 

factorized into either single- or two-spin evolution oper- 

(2) 

ators. For example, under the effect of exp {—iH^^^T) = 



N /2 (2,') 

®fe=i exp [-«/iod(j(s2fc_i,S2fe)r] (with N the number of 
spins), each pair evolves in the same way, and indepen- 
dently of all the others. 

In order to simulate the dynamics of the spin chain, we 
encode its odd- and even-numbered spins into the state 
of two physically distinguishable kinds of spin qubits, 
A and B, in the quantum hardware: S2fc-i — >■ and 
S2fc — (Fig. 1(a)). Note that hereafter we use capital 
letters to indicate the spins and times of the quantum 
simulator. 



In the latter, each pair S^S^ (S^S^^^) is 
physically connected through a spin cluster Mab {Mb a), 
whose state can be manipulated so as to effectively switch 
the coupling between the qubits. Being Mab spectrally 
distinguishable from Mb a, it will be possible to selec- 
tively switch the AB or BA couplings, still by means 

of spatially homogeneous em pulses (see below). The 

(2) 

time-evolution operator exp {—iH^^^r), is implemented 
by performing sequences of identical single- and two- 
qubit operations on the pairs S^S^. Such pulse sequence 
induces on the pair of qubits the same time-evolution in- 
duced in the time r by h^J^ on the spin pair S2fe_iS2fe. 
Analogously, sequences of pulses applied to the qubit 
pairs S^S"^j_^ reproduce the time-evolution induced on 
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FIG. 1: Top: chain of the A and B spin qubits, with inter- 
posed antiferromagnetic dimers Mab and AIba- The level 
schemes of the noninteracting units is also illustrated. Bot- 
tom: Magnetic field dependence of the energy levels for a 
A - M - B system (Eq. ©), with = 1.8, gf = 2.0, 
gf = 2.3, J = 1 cm"\ A = 30 cm"^ When the dimer 
is in the singlet ground state (four lowest levels, black), the 
qubits behave as if they were non-interacting. If the dimer is 
in triplet states (twelve highest levels, red), the energy levels 
depend on the dimer-qubit couplings. 
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FIG. 2: Simulation of the time evolution operators 
exp {—iK,siaS2p) (ctP = zz, yy, xy) for = 7r/2. Tf = 180 
ns is the duration of the pulse sequence; R, defined in the 
text, quantifies the deviation of the implemented transfor- 
mation from the target evolution operator. The em pulses 
are gaussian and linearly polarized, with a peak amplitude of 
50 G; the static field is -B = 5 T. The qubits are assumed 
to be magnetically isotropic, whereas for the dimer we set 
Qx — gS = 9y ~ Qy ~ ^- ^ is also shown for a single Trotter- 
ization step of the TIM Hamiltonian for two spins (Eq. (4)). 
Inset: time oscillations of the longitudinal average magneti- 
zation si) for the TIM with N = 3 and X = 2b. 



the spin pairs 82^824:4.1 by the terms /icvcn- The contri- 
bution of the single-spin Hamiltonian, H'^^\ is simulated 
instead by single-qubit rotations performed simultane- 
ously on all the qubits. We stress that the parameters 
defining H can be easily varied in the simulation by ap- 
propriately choosing the phases of the pulse sequence (see 
discussion below Eq. (3)). Thanks to the translational 
invariance of the simulated system, the number Nop of 
operations (i.e., pulses) does not depend on the chain 
length N . In fact. Nop is proportional to nl, where I is 
the number of terms into which e"*^"^ is factorized (see 
Eq. [T|). Besides, in many models of interest H (e.g., sym- 
metric exchange between 1/2 spins), the A and B units 
need not be physically different. 

Switching the coupling between spin qubits — Turn- 
ing on and off a physical exchange interaction between 
molecular nanomagnets on a ns timescale is presently 
unfeasible. Here we show that the effect of a switch- 
able coupling between two nearby spin qubits A and 
B can be obtained by manipulating the interconnect- 
ing spin cluster AI between them. As an illustrative 
example, we consider the case where M is a dimer, con- 
sisting of two antiferromagnetically-coupled spins (with 
Sf = Sf = 1/2, see Fig. 1). The simplest qubit- 
dimer-qubit unit is described by the Hamiltonian: 

Hamb — [Ha + Hb + Hm] + [Ham + Hbm] 

= [mbSo {gfS^ + gfS^ + g^\St^ + S^^ )) 

+ ASf .Sf] + [J5] So-Sf], (2) 

a=A,B i=1.2 



where g-^ are gyromagnetic factors, and the last term is 
the qubit-dimer coupling, whose energy scale J is typ- 
ically much smaller than /is-Bo a-nd A. A similar ef- 
fective Hamiltonian results, for ins1;ance, from two CryNi 
rings linked through a Cu2 dimer where the rings play 
the role of the spin qubits. Many generalizations of Eq. 
(2) , including less symmetric patterns of qubit-dimer ex- 
change couplings, do not alter the validity of the proposed 
scheme. The field-dependence of the energies resulting 
from Eq. ^ is depicted in Fig. 1. The four lowest 
states match those of the two isolated A and B qubits. 
In fact, as far as the dimer is in its singlet ground-state, 
{Ham) = (Hbm) — for any two-qubit state; therefore, 
A and B are effectively uncoupled, whereas they do com- 
municate if the dimer is sent by an em pulse to an excited 
state. Indeed, the twelve upper states have energies with 
fine splittings determined by J. 

This level scheme can be exploited to simulate the dy- 
namics of two generic spins induced by a Hamiltonian 
/i*-^' = KS101S213 for any choice oi a,/3 — x,y, z. The op- 
erator exp {—ih^-^^T) can be decomposed as follows: 

exp {-iKSiaS2pT) = [uia ® "2/?] e~'^^ [Ula ® U2^]^ (3) 

where A = ksizS2z, Ux = (2sy)^/^, Uy = {2sx)~^^^, and 
Uz = I. In the physical hardware described by Eq. (2), 
the single-qubit rotations Ua can be implemented by em 
pulses with frequencies Bg^ and Bgf , respectively (while 
M is left in its ground state). The two-qubit operator 
e"'^"^ would in principle require a direct interaction be- 
tween A and B. Here, instead, it's implemented by in- 
ducing in M an excitation conditioned to the state of 
the spin qubits A and B. This is obtained through two 
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simultaneous tt pulses, resonant with the gaps indicated 
by arrows in Fig. 1, followed by a repetition of the two 
pulses that bring the dimer back to its singlet ground 
state. The value of kt of the target evolution is con- 
trolled by the phase difference between the first and the 
second pair of pulses. Hence, the parameters defining 
the Hamiltonian we want to simulate can be easily var- 
ied since these merely determine the phases of the pulse 
sequences. 

In Fig. 2 we demonstrate the validity of this scheme 
by simulating the evolution operator for some represen- 
tative choices of a/3 in h^'^\ We start from the 4 pos- 
sible two-qubit basis states Itpi) (corresponding to the 
4 lowest eigenstates of Hamb, with i ~ 1,...,4), and 
we calculate their time evolution \^i(T)) induced by 
Hamb (Eq. ([2])) and by the pulse sequence. The ma- 
trix elements of the resulting transformation, Uji{T) = 
{ipj\ipi{T)), are compared with those of U{t) between 
the corresponding states of the system to be simulated: 
the distance between the two is assessed by R{T) — 
maxij \Uji(T) — Uji{T)\'^. At the end of the pulse se- 
quence (T — Tf), Uji{Tf) coincides with Uji{T), showing 
that the two qubits have actually undergone the desired 
unitary transformation (Fig. 2)'14!]. For the chosen, re- 
alistic parametrization, the duration of these simulations 
is of the order of 10^ ns. This is much less than the ex- 
pected decoherence times for optimally engineered molec- 
ular qubits, of the order of several microseconds 5]. At 
low temperatures, the coherence of each nanomagnet is 
limited by the hyperfine coupling to the nuclear spins, 
and the noise has a local character (i.e. each molecule 
interacts mainly with its own bath of nuclear spins) [isj . 
The effects of such coupling can be partially cancelled by 
spin-echo sequences. These imply the use of additional 
pulses, that can however be applied in parallel to the 
whole array, and thus independently on the system size. 
Quantum simulation of a spin- 1/2 chain — A simple 
proof-of-principle experiment can be performed by sim- 
ulating the time evolution of the transverse-field Ising 
model (TIM): 

N-l N 

'Htiu = H-nu = A ^ SkzS(k+i)z + Sfc^, (4) 
fc=i fc=i 

with Sfc — 1/2. The mapping of H^im onto the qubit 
chain is straightforward: S2k-i and S2k . 

Performing a single Trotter step of the TIM entails sim- 
ulating e~'^^'^ with K — X, followed by a rotation of two 
qubits around x by an angle br. The convergence of the 
simulated TIM evolution to the target evolution is shown 
in Fig. 2. For a generic value of A/6 the TIM evolution 
brings the system from a factorized initial state to multi- 
partite entangled ones. This is shown in the inset of Fig. 
2 for the case = 3. Here \ip{t)) starts from a ferro- 
magnetic state and it evolves passing through GHZ-like 
states [l^. This evolution is witnessed by oscillations of 
the magnetization, whose frequencies are set by the en- 
ergy gaps of the TIM. 




FIG. 3: From top to bottom: mapping of the spin-1 Hamil- 
tonian Hsi (Eq. ((5])) onto a spin- 1/2 one Hai, and encoding 
of Hsi into the spin-qubit chain ABAB .... 



Simulating S* = 1 spins and the Hubbard model — The 
simulation of Hamiltonians involving spins s > 1/2 or 
fermions requires a suitable mapping onto the qubits. 
For instance, we consider a chain of spins one (t^) with 
NN exchange interactions and single-spin crystal-field 
anisotropy: 

N-l N 

Hsi = A ^ tfe • tfc+i +Y,[d + e{tl^ - tly)], (5) 

fc=l k=l 

which reduces to the paradigmatic Haldane model for 
d = e = 0. Hsi can be mapped onto a Hamiltonian 
Hsi of 2A^ spins 1/2, with NN interactions. Indeed, 
the dynamics of a spin-1 chain is equivalent to that of 
a dimerized spin- 1/2 chain with twice the number of 
spins, provided the isotropic exchange constant (/iso) (see 
Fig. 131) is ferromagnetic and dominant. The three states 
mfc = 0,±1) of each spin tk = 1 are mapped onto 
the three triplet states of the pair of spins 1/2 (sfc,lfc), 
having total spin one. By exploiting the Wigner-Eckart 
theorem, the crystal-field terms are mapped onto axial 
(^ax) and rhombic (/rh) exchange terms: ^.^^^(sfc,!^) — 
/isoSfe • Ifc -I- laxSkzhz + Irh {skxhx " Skyhy)- The exchange 
interaction between and t^+i in Eq. ([S]) is mapped in- 
stead onto a Heisenberg coupling / between 1^ and s^+i 

(Fig. El): ft.oven(lfe,Sfc+i) = Ilk ■ Sfc+1. 

Having mapped Jigi onto a chain of spins 1/2 with 
NN interactions, we can now simulate its dynamics along 
the lines traced above. Each spin-1/2 (1^) is encoded 

(2) 

into qubit A {B) and the operators exp (— *^odd^) ^^'^ 

(2) 

exp (— i/icvcnT) are mimicked as outlined above. A simple 
proof-of-principle experiment would be the simulation of 
a single s = 1 spin experiencing tunneling of the magne- 
tization (|e| ^ |d|). For instance. Fig. |3] shows the exact 
and simulated evolution of the magnetization in the case 
d/e — 12, which can be monitored by measuring the total 
magnetization of the A — M — B system. 
The mapping of a fermionic Hamiltonian onto a 
spin one is generally nontrivial [l7| . Hereafter, we 
use the Jordan- Wigner representation in order to 
map the one-dimensional Hubbard model, Tiuuh = 
-^H Y.kcri4:a'^k+ia + h.c.) + U Y,k ^k^'^'ki, onto a chain 
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FIG. 4: Line: exact time evolution of {t^) for a single t = 1 
spin with d/e — 12, Eq. ([5]). The state oscillates between 
\m = 1) and \m — —1) due to quantum tunneling across 
the anisotropy barrier. Points: time evolution simulated by a 
single A — M — B unit initially prepared in its ground state. 
The pulse sequence is set by the mapping of Eq. ((5]) onto 
the Hamiltonian of two 1/2 spins (Fig. 3). We plot the z- 
component of total spin of the A — M — B unit, which could be 
easily extracted by measuring the magnetization of a crystal 
of noninteracting units. The duration Tf of the pulse sequence 
implementing the simulation is about 480 ns, independently of 
the simulated time t. Note that to perform each simulation 
(i.e., to extract each point) the A — AI ~ B unit has to be 
reinitialized to its ground state. 

of 1/2 spins 0: 

N-l 

fe=l a=x,y 
N N 

+ U"^Skzlkz + U/2j2iskz+lkz) (6) 

fc=l k=l 

where Sk and are two families of spin-1/2 operators. 
Unlike the previously considered cases, i?Hub is not a 
one-dimensional Hamiltonian with NN interactions only, 
as each spin couples to three other ones. This requires to 
proceed in two steps: initially we encode into the qubit 
pairs (S^,Sf) the spins (lfe,Sfe) for even k and {sk,lk) 
for odd k , respectively (see Fig. [SJ. The couplings be- 
tween pairs of spins that are encoded into neighboring 
qubits can be simulated as for the previously considered 
models. These couplings include the two z terms in Eq. 

and half of the xy terms, i.e., the even transverse 
s — s bonds and the odd transverse I — I bonds. In or- 
der to simulate the remaining transverse two-spin terms 
in i?Hubj corresponding to couplings between qubits that 
are initially third nearest-neighbours, we swap the state 
of all the S^Sf pairs. This SWAP gate can be performed 
by the same method used to implement two-qubit opera- 
tions. Third nearest neighbours in -ffnub now correspond 
to NN in the simulator, and the evolution induced by the 
associated couplings can be simulated exactly as above. 

The feasibility of our scheme with available technol- 
ogy relies on the lack of local-control requirements, as 
only uniform em pulses are involved. We have illustrated 
the simplest possible implementation of the idea, work- 
ing for uniform or ^ — _B 1-dimensional Hamiltonians. 
However, extensions of this approach allow the simula- 
tion of a much larger class of Hamiltonians, including 
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FIG. 5: Mapping of the Hubbard model Huuh onto the spin 
Hamiltonian -ffnub involving the Sk and Ik 1/2 spins. 



higher-dimensional ones. One possibility is to use spin- 
qubit arrays that reproduce the dimensionality D of the 
system to simulate, with interposed nanomagnets M to 
switch the interaction between adjacent qubits. For in- 
stance, in a square lattice U{t — >■ 0) is first decomposed 
by the Trotter formula into two evolution operators de- 
scribing a collection of identical chains (along x and y). 
The couplings within the chains can then be simulated in 
parallel by the method described above by using an ar- 
ray of nanomagnets with rectangular symmetry. In fact, 
in order to selectively address the x and y chains, the 
nanomagnets M switching the interaction along one di- 
rection need to be spectrally distinguishable from those 
operating on another direction. 

Alternatively, it is possible to develop other simulation 
schemes keeping a one-dimensional topology of the hard- 
ware and the use of uniform pulses, at the cost of a more 
complex and less parallel algorithm. The limiting case 
(where all two-body terms in the Trotter-decomposed 
time evolution are implemented sequentially) is repre- 
sented by a scheme where half of the molecular qubits 
are used as auxiliary units instead of logical qubits. The 
Z?-dimensional target Hamiltonian is mapped onto a ID 
Hamiltonian with long-range couplings. Since the latter 
prevent adopting the scheme described above for nearest- 
neighbor ID Hamiltonians, the various two-body terms 
in the target Hamiltonian are simulated sequentially by 
the use of a control unit (in the spirit of [19]). Conditional 
excitation of the interposed dimers remains a key ingre- 
dient to induce the desired evolution, and the auxiliary 
units are exploited to attain local control with uniform 
pulses 20] . 

The capability of simulating Hamiltonians can also be ex- 
ploited in order to map experimentally accessible quan- 
tities (O) onto a set of non-accessible observables (C = 
^iHtQ^-iHty Observables corresponding to (sums of) 
single-qubit terms can be mapped, e.g., onto pair cor- 
relation functions in the case where H corresponds to 

(2) (2) 

a dimerized system {H — H^^'^ or H — Hcvcn), or to 
higher-order correlation functions in the case of a more 
general H^. 

In conclusion, we have shown that arrays of molecular 
nanomagnets can be used as quantum simulators of dif- 
ferent model Hamiltonians with translational invariance 
and short-range interactions. We have proposed proof- 
of-principle implementations, where the means required 
for manipulating the system and measuring the relevant 
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observables can be provided by the current technology. 
Nanomagnet-dimer-nanomagnet supramolecular trimers 
and chains which can be exploited to implement our 
scheme are also currently being synthesized ^2l] . 
We thank M. AfFronte for useful discussions. We ac- 
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